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Abstract 

In finite volume the partition function of QCD witli a given is a sum of different topological 
sectors with a weight primarily determined by the topological susceptibility. If a physical observable 
is evaluated only in a fixed topological sector, the result deviates from the true expectation value by 
an amount proportional to the inverse space-time volume 'i/V . Using the saddle point expansion, 
we derive formulas to express the correction due to the fixed topological charge in terms of a 1/V 
expansion. Applying this formula, we propose a class of methods to determine the topological 
susceptibility in QCD from various correlation functions calculated in a fixed topological sector. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) in four space-time dimensions allows topologically 
non-trivial gauge configurations labeled by a winding number or a topological charge Q. 
The path integral to define the partition function of QCD includes an integral over configu- 
rations with arbitrary Q. In order to ensure the cluster decomposition property of physical 
observables, the weight among different topological sectors must be e*^*^, which defines the 9 
vacuum of QCD. If one considers a path integral restricted in a fixed topological sector, the 
cluster decomposition property — one of the necessary properties of quantum field theory — 
is violated In this paper we address the question of fixing topology in the context of 
non-perturbative calculation of QCD on the lattice. But the analysis does not depend on 
any particular regularization of the theory. 

In the lattice QCD simulations the inclusion of the effects of dynamical quarks is com- 
putationally most demanding. Since the direct computation of the fermion determinant 
det(D + m)^, with D the Dirac operator on the lattice and m the quark mass, requires 
prohibitive computational cost, one usually introduces a pseudo-fermion field (f) to write the 
determinant in the form J [c?0'''] exp(— + m)~^0), so that the problem is reduced to 
an evaluation of the inverse fermion matrix {D + m)^^. Since the effective lattice action 
becomes non-local, the Monte Carlo updation is most efficiently done by updating all the 
gauge links on the lattice at the same time introducing a molecular dynamics evolution. The 
most popular such algorithm to date is the Hybrid Monte Carlo algorithm that combines 
the molecular dynamics evolution with a Metropolis accept/reject step. 

With the molecular dynamics evolution, the ergodicity becomes a potential problem when 
there exist more than one regions of phase space that are separated by some potential wall, 
because the "kinetic energy" of the molecular dynamics system may not be enough to go 
through the potential wall. This situation happens for QCD in four dimensions because of the 
non-trivial topological sectors. In the continuum theory the potential barrier is infinite and 
the gauge configurations in different topological sectors cannot be reached by a single stream 
of the continuous evolution. On the lattice, the potential barrier is finite (order of inverse 
lattice spacing 1/a) and the probability of tunneling among different topological sectors is 
non-zero, but will exponentially drop (~ e~^°/" with U^/a a, nominal potential height) as 
the continuum limit is approached. This means that the correct sampling of topological 
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charge and thus the vahd simulation of the ^-vacuum of QCD will become increasingly more 



difficult 



a 



In fact, this problem already manifests itse 



fermion simulations by the JLQCD collaboration 
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f_in the recent dynamical overlap 
since they explicitly introduce 



a term that prevents topology change ^ in order to avoid large numerical cost due to the 
discontinuity of the overlap operator along the topology boundary. 

One may then ask whether lattice QCD simulations with a fixed topological charge are 
useful in general to extract physics of real world, i.e. QCD at a given value of 9. (If the CP 
symmetry is preserved, ^ = 0. Small but non-zero 6 implies an interesting physics case to 
give rise to the neutron electric dipole moment.) The fixed topology simulation is desirable 
in the study of the e-regime of QCD, since the physical quantities have striking dependence 
on the topological charge, which is an important part of the physics we are interested in. 
(For a recent unquenched lattice simulation in the e-regime, see 10|, lUl-) But in the p- 



regime, the fixed topology simulation gives rise to a systematic effect. In this paper we try 
to answer this question theoretically without relying on any particular model or assumption. 
Extending a previous work by Brower et al. 12], we show that observables calculated at a 
fixed topological charge are different from those at ^ = by an amount of order of inverse 
space-time volume l/V , so that both agree in the infinite volume limit. 

This result can be understood intuitively within the instanton picture of the topological 
excitation in the QCD vacuum. The instanton is a local object carrying a unit topological 
charge. Starting from a trivial topological sector, an instanton and anti-instanton can pair- 
create or -annihilate without changing the net topology and then may apart from each 
other. If we look at a region (or sub- volume) including the (anti-) instanton, that region 
has locally non-trivial topological charge. In this way, the topological fiuctuation may occur 
even if the net topological charge is kept fixed. A supporting evidence of this expectation has 
been found in the chirality density of low-lying eigenmodes of the overlap-Dirac operator, i. e. 
local chirality is strongly peaked at ±1 even on the gauge configurations with net topological 
charge Q = . The effect of fixing the net topology will become relatively small as the 
total volume becomes larger, as there is more chance to create (or annihilate) the instanton- 
anti-instanton pairs. The topological susceptibility Xt characterizes the magnitude of these 
local topological fiuctuations. 

In fact, in our formulation the topological susceptibility enters in the difference of 0{1/V) 
between the fixed topological sector and the fixed 9 vacuum. This type of finite volume effect 



3 



can even be evaluated, once we know the coefficient of l/l^, which is found in general to 
be a second derivative of the physical quantity of interest with respect to 6*, as well as the 
topological susceptibility Xt- In this paper, we describe the method to extract Xt from a 
gauge ensemble in the fixed topological sector, leaving the estimate of the finite size effects 
for other physical quantities for a future publication. 

This paper is organized as follows. In Section [Tll by using the saddle point expansion, 
we derive a formula which express the partition function at a fixed topological charge in 
terms of the partition function at 6' = and its derivatives with respect to 6. We extend 
this formula to the case for arbitrary correlation functions. From these results it is easy 
to show that arbitrary correlation functions at fixed topological charge agree with those at 
^ = in the infinite volume limit. In Section [Till as an application of the above formulas, 
we show that the topological susceptibility Xt can be extracted from two-point function of 
the topological charge density at a fixed topological charge. Through the Ward-Takahashi 
identities, we then relate the two-point function of the topological charge density with a two- 
point function of the pseudo-scalar density, which is more suitable for actual calculations 
with the overlap fermions. An extension to the case of four- and three-point correlation 
functions is also discussed. In Section llVt we consider a case of CP-odd observables taking 
the calculation of the neutron electric dipole moment as an example. Our conclusion is 
given in Section |Vl In Appendix |X] we check a validity of the saddle point expansion, by 
comparing the expansion with an exact calculation for a simple example. Results at order 
are summarized in Appendix [Bl 

II. GENERAL FORMULA 

A. Partition function at a fixed topological charge 

We ffist consider the partition function in the 6 vacuum defined by 

z{e) ^ {e\e) = exp[-vE{e)], (i) 

where E{6) is the energy per (space-time) volume V. Without loss of generality we take 
Z[0) = 1 as a normalization, which is equivalent to E{0) = 0. E{6) satisfies the conditions 
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E{9 + 2ti) = E{9) = E{—9). The topological susceptibility at ^ = is defined by 

X, - (2) 

which is obtained by a second derivative of E{9) 

£E{9) 



Xt 



d9^ 



(3) 

Since ^ by definition, ^ = is a local minimum of E{9). Moreover, Vafa and Witten 



proved that Z{0) > Z{9) [l3|, which leads to E{0) < E{9) for ^ 0. Namely, ^ = is 
the global minimum of the function E{9). While we assume that E{9) is analytic at 6' = 0, 
E{9) may have a non-analyticity at ^ 7^ 0, in particular at ^ = tt. For instance, the Chiral 
Perturbation Theory (ChPT) at the leading order gives 

E{9) = XtN^f[l-cos{9/Nf)l -Tx<9<7i (4) 

for Nf the number ^ flavors. A 9 derivative of this particular form is discontinuous at 



9 = ±7r for A'^^ > 1 {l^. From these general conditions we can expand E{9) around 9 = 0, 

00 

n=l ^ 

The partition function at a fixed topological charge Q is a Fourier coefficient of the 
periodic function Z{9) 

Zq = ^ rd9Z{9)exp{z9Q) = ^ F d9 exp{-VE{9)), (6) 

where F{9) = E{9) — i9Q/V , and the Fourier expansion of Z{9) is written as 

Q=0,±l,±2,- 

For a large enough volume, we can evaluate the 9 integral in (jS]) by the saddle point expan- 
sion. The saddle point 9^ is given by 

9. = ^^{l + 0{5')) (8) 

where 5 = Q/ixtV)- We then expand F{9) as 

m = Fi9^) + -^m9 - 9^)' + ^ -^{e - 9^, (9) 

n=3 



where -E^"^ is the n-th derivative of E{9) with respect to 9 at 9 = 9^, and is given by 



VF{9,) 
E^'\9,) 



9,C2n{l + OiS')). 



(10) 

(11) 

(12) 
(13) 



Hereafter we omit an argument 9c of E^"'\9c) and simply write it as E^"^^ unless otherwise 



stated. By a change of variable s = V E^^'>V{9 — 9c) we can rewrite the integral as 



(is exp 



'(2)y{-7r-9c) 



71=3 



n! 



(is exp 



S2 ^ 



n=3 



V27rE(2)1/ 



exp 



E 

n=3 



n! 



hO(e-). 



We defined 



-L / dse-'"'''f{s) 
v27r j_oo 



(14) 



(15) 



with which {s^^) = {2n — 1)!!. Neglecting exponentially suppressed terms and expanding in 
powers of 1/V, we obtain 



V27rEW 
1 

, exp 



E 



1 n=3 
B(E(2)y)2 



2xtV^, 



1 - 



C4 



y2' 



5^ 



(16) 



This shows that, as long as 5 ^ 1 (equivalently, Q -C Xf^)? the distribution of Q becomes 
the Gaussian distribution. Note however that the distribution can deviate from a Gaussian 
for Q = 0{V), as numerically observed on a quenched lattice 
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B. Correlation functions 



Here we consider an arbitrary correlation function 

G{9) = {9\Oi02---0^\9), 
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(17) 



whose Fourier coefficient at a fixed topological charge Q is defined by 



Gq = ^^ / dez{e)G{e)eM^eQ). 



(18) 



Note that the operators Oi do not contain 9, and the 9 dependence comes solely from the 
vacuum angle. Using the saddle point expansion as before, we obtain 

X e-VF{0c) /-oo 



Gc 



ds exp 



' 9 2-^ 



2 ^ n! y^WW 



G 9 



(19) 



where we define 



iifm 



I ^-VF{e,) /-oo 



ds f{s) exp 



s 



E 

n=3 



(20) 



and G^"'^ is the n-th derivative of G with respect to 9. Using the formulas 



((1)) 
iis')) 

{{s)) 



2(^(2)1^)2 4(E{2)V^)3 

3!! + 0(^-1), 

~2(E(2)\/)3/2 



+ o(\/-2), 

5(^(3)^)2 



(21) 
(22) 
(23) 



+ 



3(E(2)\/)2 4(E(2)y)3y 8(E(2)\/)5/2 



+ 0(^-^/2) ^24) 



■2(E(2)\/)3/2 



+ o(\/- 



-3/2N 



(25) 



we finally obtain 

Gq = Gi9,) + G^^\9,) 

-G^'\9c 



2E(2)y 



E^^'^V ^ 5(E(3)\/)2 



2(E(2)V^)^ 



2(E(2)F)2 4(E(2)\/)3 

4E(4)y 5(^(3)^)2 



-G^^\9, 

3(E(2)y)2 ' 4(E(2)V)3; ' 8(E(2)r) 



-G^^H^c)jS^ + 0(V-) 



12(E(2)y)3 

The above expansion is valid for any 9c as long as G^^\9c) /G{9c) = 0(1) as ^ oo. 
Depending on the size of 9^, we can further expand the above formula. If we take 9, 



(26) 
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0(y (equivalently, Q = 0(1)), we have 
Go 



G(0) + GW(0)^, + G(2)(0)| + G(^^(0)^^ 



+GW(0)- 



1 



G'(0)+G(')(0) 



(E(2)y)2 

1 



G(^)(0) 
1 



2(E(2)r)^ 



3(^(2)^)27 8(^(2)^)- 



C4 



+0^(0) 



1 



C4 



G^^^HO)- 



+ G(^)(0) 



1 



If G is CP-even, G is an even function of 9, so that 



^even 
Q 



G(0) + G(2)(0) 



1 



1 - 



C4 







0{V-^). 



(27) 



while, if G is CP-odd, we have 



^odd 



G^'\0) 



C4 



G(=^)(0)- 



(28) 



(29) 



In order to claim that the above expansion is converg ent, G(")(0)/G(^) must be 0(1). This 
condition is satisfied because the expansion 



G'(^) =G'(0) + ^G(")(0) — 



(30) 



n=l 



is vahd for = 0(1). 

The formula ( l28l) provides an estimate of the finite size effect due to the fixed topological 
charge. The leading correction is of order 0(1/^) as advertised. The dimension is compen- 



sated by the topological susceptibility Xt- In ChPT it is evaluated as Xt = mTs/Nr 



3y the 



sea quark mass m and the chiral condensate S as well as the number of flavors Nf [19|] . The 
flnite volume correction is suppressed when the quark mass is larger than 1/(SV), while 
becomes signiflcant when m ~ l/(Sy). This is consistent with the fact that the topological 
charge has a strong effect in the e-regime, which is characterized by mT,V < 1. The cor- 
rection has a coefficient 0(^^(0) that represents the ^-dependence of the correlator. This is 
not known in general but can be fitted with lattice data at various Q. In the mass region 
where ChPT is applicable, it can also be estimated, as done in [l^ at the tree level for 
pseudo-scalar meson mass. One-loop calculations are in progress. 
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The other interesting formula (129|) suggests a possibihty to calculate CP-odd observables, 
such as the neutron electric dipole moment, at a fixed non-zero topological charge, which 
will be discussed in Section [IVl 



III. TOPOLOGICAL SUSCEPTIBILITY 

In this section, we propose methods to extract the topological susceptibility Xt from 
correlation functions at a fixed topological charge Q. 



A. Two-point correlation function 

1. Bosonic formula 

Suppose that there is a well-defined local operator uj{x) that measures the local topological 
charge. The global topological charge Q is then obtained as Q = / d'^xuj{x), and the 
topological susceptibility is Xt = J d'^x{uj{x)uj{0)), where the expectation value is taken for 
the 6* = vacuum. 

Since uj{x)uj{0) is CP-even, Eq. ( l28l) gives 

{uj{x)uj{0))q = {uj{x)uj{0)) + (^(x)a;(0))(^); ^ 



2Vxt V Vxt 2xiV 
+ (.;(xV(0))(^)^ + O(V-^). (31) 

Here (O) means a vacuum expectation value of (9 at 6' = and (O)*^") denotes its n-th 
derivative with respect to 9. In the large separation limit |x| oo, the CP invariance at 
6* = and the clustering property at a fixed 9 [l| gives 

(cj(x)cj(O)) ^ (cj(x))(cj(0)) = 0. (32) 

In addition, 

(^(x)^(0))(2) = -(^(x)^(0)g2) + (^(a;)^(0))(g2), (33) 

whose second term vanishes as |x| oo. Denoting a connected VEV as {■ ■ ■)ci the first 
term can be written as 

^ d'^y d'^z{uj{x)uj{Qi)uj{y)uj{z)) = jd^yd'^z {uj{x)uj{Q)uj{y)uj{z))c 

+ {u:{x)ujm{uj{y)uj{z)) + {uj{x)u{y)){uj{Q)uj{z)) + {uj{x)u{z)){uj{y)ujm 
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As |x| — > oo, the first and the second terms vanish and we obtain 

U% jd^z{uj{x)uj{y)){uj{z)uj{'d)) + {y^ z) = 2 Jd^y{iu{y)uj{0)) jd'^z{uj{x)uj{z)) . (34) 

Using the translational invariance, this becomes 2%^. Therefore, the term (co'(x)co'(O))*-^'' in 
f lHT]) gives — in the large separation limit. 
Similarly, we consider 

(a;(x)a;(0))(^) = {uj{x)uj{Q)Q^) - Q{uj{x)uj{{))Q^){Q^) + (cu(x)a;(0)){6(g2)2 _ (g4^|^ (35) 
The first term is written as 

{uj{x)uii{{^)Q'^) = I '^^d^Xi{uj{x)uj{Q)uLi{xi)uj{x2)uj{x^)uLi{x^)^ (36) 

i 

which may be decomposed in terms of the connected VEVs. Since the terms containing 
both uj{x) and uj{0) in the same connected VEV vanish, only the following terms remain: 

Y\ d'^ Xi %{u}{x)u){xi)u}{x2)oj{x^))c{oj{Q)uj{xi)) 



+l2{u{x)u{xi)){u{x2)uj{x^)){u{Q)uj{x^))\ = 8xt^ + l2Vxl (37) 
Putting this into ( l35l) we have 

{u{x)um^'^ - 8x.^ = ^^log^W = -8X*C4 (38) 

for the fourth-derivative term in (13T1) . 
Gathering all terms we arrive at 

hm (c.(x)c.(0))q = yA%-^^- ^1 + ^(^"') + 0(e-"VN), (39) 

I^Hlargc V \ V 2XtV J 

where the flavor singlet pseudo-scalar meson mass, m^/, is the lightest mass of possible 
intermediate states. This result explicitly shows that QCD at a fixed topological charge 
(in a finite volume) is sick as a quantum field theory, because the clustering property is 
violated. However, the magnitude of the violation can be estimated as a 1/V expansion 
without introducing strong assumptions on the details of the dynamics of QCD. The same 
result (up to the C4 term) was obtained in a context of a two-dimensional model in [2^ . 

Physical quantities such as the topological susceptibility xt can be obtained through 
( l39l) . In practice, this formula will be used for a finite separation x instead of |x| 00. 
The clustering property (l32l) in the 6 vacuum receives a correction of order of e""*''''^', which 
vanishes quickly because the flavor singlet meson rj' acquires a large mass due to the axial 
anomaly of QCD. 
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2. Fermionic formula 



We now express the bosonic correlation function {uj{x)uj{0)) in terms of a fermionic one 
using the anomalous axial U(l) Ward-Takahashi (WT) identities for an arbitrary operator 
O: 

{df,A^{x)0 - 2mP{x)0 + 2uj{x)0 + 6^0) = 0, (40) 

where Afj_{x) = X]/- ^'^(^)7m75V^'^(^) P{x) = J2 r '^'^ i^) ^^e flavor sin- 
glet axial-vector current and pseudo-scalar density, respectively, and 5^0 denotes a axial 
rotation of the operator O at x. Here the quark field has a flavor index / running from 1 
to Nf. The expectation value (■ ■ ■ ) in ( l40l) can be taken either in the 9 vacuum or in the 
fixed Q sector, since the WT identities are valid with any external states. Combining the 
following two WT identities (for O = 2mP(0) and O = 2uj{x) with x ^ 0) 

(2mP(x)2mP(0)) = {d^A^{x)2mP{0)) + (2cj(a;)2mP(0)) (41) 
{2u{x)2m.P{0)) = {2uj{x)d^A^,{0)) + {2u{x)2uj{0)), (42) 

we obtain a relation at a fixed topological charge 

(2mP(x)2mP(0))Q = (2w(x)2a;(0))Q + {d^A^{x)2mPiO))Q + {2u;{x)d^A^{0))Q. (43) 
In the large separation limit |x| — * oo, using (l28l) the second term becomes 
{d^A^{x)2mP{0))Q ^ {d^A^ix)){2mP{0)) 
+ 2!7^ - C " ^) + ^2(9,M^)^rnP{0)r. (44) 

The first term on the right hand side vanishes, and the second term is evaluated as 

(a^A^(x)2mP(0))(') = {d^A^ix)2mP{0)Q^) ^ 2{d^A^ix)Q){2mP{0)Q). (45) 
Since the translational invariance (plus an appropriate boundary condition) leads to 

{d,A^{x)Q) = d^x{d^A^{x)Q) = (46) 

we obtain {df^A^(x)2mP{0))^^^ = in the large separation limit. 
We next consider the fourth derivative term. In general, 

(A(a;)P(0))(^) = {A{x)B{0)Q') - 6{A{x)B{0)Q'){Q') + {A{x)B{0)){6{QY - {Q')} (47) 
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when {A{x)) = (-B(O)) = 0. Using relations 

{A{x)B{0)Q^) ^ 2{A{x)Q),{B{0)Q),, (48) 
{A{x)B{0)Q^) = {A{x)B{0)Q'),+ {A{x)B{0)),{Q'), + 6{A{x)B{0)Q'),{Q'), 

+3{A{x)BiO)),{Q')l + 4{A{x)Q)c{B{0)Q'), + 4{B{0)Q)M{^)Q')c 

+12{A{x)Q),{BiO)Q),{Q''), 
^ A{A{x)Q),{B{0)Q^), + 4{B{0)Q)M{x)Q^)c 

+ 12(A(x)g)e(S(0)g)e(Q')e, (49) 

the general formula reads 

(A(x)5(0))(^) ^ 4(A(x)Q),(fi(0)g=^), + 4{B{0)Q)Mi^)Q^)c. (50) 

For yl(x) = dfj,Afj_{x), the translational invariance again gives {A{x)Q)c = {A{x)Q^)c = 0, 
and therefore (9^A^(x)2mP(0))*^'^) — 0. We thus conclude that the term (9^A^(x)2mP(0))Q 
on the right hand side of (1151) vanishes in the large separation limit |x| oo, and the same 
conclusion also holds for the last term {2uj{x)df^Afj_{0)) q in (1431) . 
Combining these we finally obtain 

hm (mP(x)mP(0))Q = lim {uj{xMO))q = U^-Xt- ^) + O(e--vl-l). 

>large >largc V \ V '^Xt^ J 

(51) 

Namely, the topological susceptibility can also be calculated with the singlet pseudo-scalar 
density correlator at a fixed topology. If we apply the above evaluation to the non-singlet 
WT identities, we can show 

lim (2mP'^(x)2mP''(0))Q = lim {d^Al{x)2mP\<d))Q = (52) 

|a;|— >largc |a:|— >large 

where A'^^{x) = tlj{x)'~ffj,'~f5T"-ip{x) and P°' = il!{x)'-f^T°'ilj{x) are the flavor non-singlet axial- 
vector current and pseudo-scalar density, respectively, with tr = and tr(T"T^) = 6"'''/Nf. 
By subtracting this equation from (15T]) . we may write 

lim (mP(x)mP(0))f ^ = lim [{mP{x)mP{0))Q - (mP''(x)mP''(0))Q] 

>large |a::|— ♦large 



where the superscript "disc" represents a disconnected contribution, and no summation 
is taken here for a. Since the large separation limit of the non-singlet correlator fl52|) is 
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saturated only after a slowly dumping factor e '"'"'^l with the pion mass m.^, the singlet 
correlator (l5Ti) is preferable for a practical use. 



B. Four-point correlation function 

Here we extend our analysis to a four-point function. 

1. Bosonic formula 

We consider a CP-even observable 

{u{xi)uj{x2)uj{x3)u{x4))q = {uj{xi)uj{x2)uj{x3,)uj{x4)) 



+ (CJ (Xi ) (X2 ) (X3 ) (X4 ) ) — ^ 



2Vxt V 2xiv 

+ {uj{x,Mx2P{xsP{x,)y'^-^ + 0(1^-3) (54) 

We assume that four topological charge density operators are far apart from any oth- 
ers {\xi — Xj\ large for Vz Vj). In this limit both {Lj{xi)uj{x2)uj{x3)Lj{xi)) and 
{uj{xi)uj{x2)uj{x3)uj{x4))^'^^ vauish and {uj{xi)uj{x2)io{x3)uj{x4))^^^ becomes 24%^, taking the 
procedure used to obtain the analogous formula (138|) in the two-point function case. We 
obtain 



2 



V 

{u{xi)uj{x2)uj{x3)u{Xi))Q = 



1 



-I 2 



1 + ^(C4-Q Xt 



OiV-^). (55) 



In addition to the leading 0(1/^^) term, we also give a next-to-leading 0{1/V^) correction. 
Details of this 0{1/V^) calculation are given in Appendix [Bl 

2. Fermionic formula 

Again we use the anomalous axial WT identities, which give a chain of equations 

{P1P2P3P,) = {dA,P2P3P,) + {001P2P3P,) (56) 

{UJ1P2P3P4) = {dA2UJlP3P4) + {tOlUJ2P3PA) (57) 
{(jJlUJ2P3Pi) = {dA3UJiUJ2Pi) + {uJiUJ2U}3Pi) (58) 
{UiUJ2UJ3Pi) = {dA4UJiUJ2UJ3) + {uJiUJ2^^3^4:) (59) 
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where we have used a short-hand notation, Pi = mP{xi), dAi = d^A^{xi) /2, and uoi = uj{xi). 
The operators are far apart from others as before. Since the WT identities are vahd in the 
fixed topological sector, we obtain 

(PiP2^3A)q = {uJiUJ2UJ3UJ4)q 

+ {dAiP2PzPi + uJidAiP^Pi + u}iLJ2dA^Pi + LJiLJ2i^^dAi)Q. (60) 

From the clustering property and CP invariance the correlator (dAiOjki) vanishes, where 
Ojki = FjF^Fi with Fi = Pi or Ui [i, j, k, I are all different). In addition, 
{dAiOjki)^"^^ vanishes, because it is written as —{dAiOjkiQ^) — {dAiQ) {O jkiQ) ■ Simi- 
larly, (dAiOjki)^^^ vanishes, because it is {dAiOjkiQ^) 4:{dAiQ){0jkiQ^) and {d^A{x)Q) = 
(l/V) f d'^x {d^A{x)Q) = 0. We therefore obtain {dAiOjki)Q — > 0, which leads to 

2 

+ 0{V-^). (61) 



xtv 



Similarly the non-singlet WT identities imply that {P^P^PkPi)Q and {P^P^P^Pi)q van- 
ish in the large separation limit, where P"" = P°'{xi) is a non-singlet PS density and i, j, /c, 
/ are all different from others. We then conclude that the disconnected correlator has the 
same asymptotic form as {ijj\UJ2ijJzUJA)Q: 

{P,P2P,P,)t'' = (PlP2P3P4)Q-{(PfP2"P3P4)Q + (PlP2P3"Pr)Q-(PfP2"P3'P4')Q} 

-{2 ^ 3} - {2 ^ 4} 

0{V-'), (62) 



Q Xt 



1 



l + -3Y7(c4-Q'Xt) 



where no sum is taken for a^h. 

The formulas (l55l) and ( 1621) provide another method to extract the topological suscepti- 
bility from a fixed topological sector. Since the contamination from the higher order effect 
in 1 jV with the unknown constant C4 appears with a different coefficient from the two-point 
case (!39|) and (1511) . both xt and C4 can be extracted by combining the two- and four-point 
correlators. Note that this calculation is free from short- distance singularities, because the 
operators are explicitly put apart from others. This is in contrast to the usual definition 
Xt = (l/V) J c/^x(cu(x)cu(0)) (at a fixed 6). 
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C. Three-point correlation function 



Topological susceptibility can be extracted also from CP-odd observables. The simplest 
case is given by a three-point function as 

iQ 

lim {uj{xi)uj{x2Mxs))q = lim {uj{xiMx2)uj{x3)y^^ + 0{V-^) 

Fi-a;j|— >large \xi-Xj\^la.vgc Z,\^V 



Xj |—> largo 



-3Xi 



Q_ 



1 + 



7C4 



0{V-^) 



(63) 



As in the cases of two- and four-point functions we may rewrite the above result in terms 
of the fermionic quantities 



lim {mP{xi)mP{x2)mP{xs))Q = lim {mP {xi)mP{x2)mP [x-^)) A 

|xi—a;j|—> large Ix^— |— >largc 



disc 



lim {u{xi)u{x2)u}{x^))c 

Xj |—> large 



Q_ 

^2 



7ca 



0{V^^). (64) 



As expected, the three-point function is useful only when the net topological charge is non- 



zero. 



IV. CP-ODD ELECTROMAGNETIC FORM FACTOR 



The three-point function considered above is an example of CP-odd observable. In our 
formulation, as the formula (129!) implies, the CP-odd observable (or its first derivative with 
respect to 6) can be calculated in a (non-zero) fixed topological sector. Another interesting 
example of this class is the neutron electric dipole moment. 

The CP-odd electromagnetic form factor F^{q^) is related to the neutron electric dipole 
moment as 



d 



N 



q^->o 2m 



N 



(65) 



where g is a momentum transfer. It has been shown in 2l| that F^i^q"^) can be extracted 
from the correlation functions {N{p)J^^{q)N{p')Q) and {N{p)N{p)Q) , where N{p) is a 
neutron interpolating field and J™ is the electromagnetic current. The momentum transfer 
is denoted hj q = p — p'. 
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Once Xt is extracted from the methods explained in the previous section, these neutron 
correlators can be extracted from at a fixed topological charge using fl29|) 

(iV(p)^''(?)iV(p'))?' = {N{p)Jr(^)N{p')Q)^ + 0{V-') (66) 
{N{p)N{p)rQ'^ = {Nip)N{p)Q)^ + 0{V~'), (67) 

where the superscript "odd" is understood that a CP-odd part of the correlation function 
is considered. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have derived general formulas which express arbitrary correlation func- 
tions at a fixed topological charge Q in terms of the same correlation function (and its 
derivatives) in the 6 vacuum. As expected from the intuitive argument of local topological 
excitations, the difference between the fixed Q vacuum and the fixed 6 vacuum disappears in 
the large volume limit as 1/V. The relation is established only using fundamental properties 
of the quantum field theory, such as the cluster decomposition principle. 

These formulas open a new possibility to calculate physical quantities in the lattice QCD 
simulations at a fixed topological charge. This will become unavoidable as the continuum 
limit is approached, irrespective of the lattice fermion formulation one employs, as far as the 
algorithm is based on the continuous evolution of the gauge field. (A proposal to avoid this 



limitation has recently been proposed [22| . Its numerical feasibility is yet to be investigated.) 

For the correlators of the local topological charge operator, the 1/V expansion is worked 
out for two-, three-, and four-point functions to the second non-trivial order. The local 
topological charge operator may either be a bosonic one or a fermionic singlet pseudo-scalar 
operator. At the leading non-trivial order the topological susceptibility Xt appears as an 
expansion parameter, and a higher order parameter C4 appears at the second order. In prin- 
ciple, these parameters can be determined by the lattice data. The different correlators and 
different V and Q may be used to check the results. This method is free from short-distance 
singularities, since the local topological charge operators are put apart from others and no 
contact term appears. Numerical calculation is in progress by the JLQCD collaboration 
on the gauge configurations generated with dynamical overlap fermion js], O, [?, [s]. Once 
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these parameters are numerically obtained, they can be used as input parameters for other 
physical observables. 

The limitation of the formulas comes from the use of the saddle point expansion. It 
requires that the volume is large enough, XtV ^ 1) that the local topological fluctuation is 
in fact active. It corresponds to the condition (0|(5^|0) ^ 1 for that volume. In addition, 
in order that the saddle point 6c — Q/ixtV) can be expanded around ^ = 0, the (fixed) 
topological charge \Q\ must be much smaller than xtV = (0|(5^|0). Since Xt is estimated 
in ChPT as Xt = rnL/Nf for Nf flavors of sea quarks, the condition is \Q\ ^ mSV . This 
means that the system must be well apart from the e-regime [rnLV < 1). This makes 
sense, because in the e-regime the physical quantities have substantial dependence on the 
topological charge and can not be simply expressed by the saddle point expansion. 

Our general formulas can also be applied to other observables, such as masses, decay 
constants, and matrix elements. Although the knowledge of the 9 dependence of the observ- 
ables is needed, only a few derivatives with respect to 9 are sufficient for small ^, which is 
physically most relevant. Using the systematic 1/V expansion one can in principle extract 
those derivatives to arbitrary finite order by looking at the Q depenence. Combining them 
with Xti C4; the physical observables in the small 9 vacuum can be reconstructed. This 
is not surprising, because if we were able to determine all the coefficients in E{9) we could 
compute observables in the 9 vacuum by a reweighting method. The point of our work 
is to present a practically feasible strategy in which the computational effort is drastically 
reduced. For the quantities well described by ChPT, estimation of the 9 dependence based 
on the chiral lagrangian is possible. Some results at the leading order were obtained in 



121, and we are extending them to the next-to-leading order. This can also be used as an 



independent consistency check. 

Once Xt is extracted, we can calculate the first derivative of the CP-odd quantities from 
simulations at a fixed non-zero topological charge using ( l29l) . The most interesting such 
quantity in the context of lattice QCD calculation is the neutron electric dipole moment. 
For this quantity, the calculation of the CP-odd form factor in the fixed non-zero topological 
charge suffices to predict the physical result in the 9 vacuum. 
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APPENDIX A: COMPARISON OF THE SADDLE POINT EXPANSION WITH 
AN EXACT e INTEGRATION 

In this appendix we demonstrate that the saddle point expansion reproduces the results 
from the exact 9 integration, by considering a simple model for E{9). 
We consider the following form of the 9 dependence: 

E{9) = xt[l-co^9]. (Al) 

The partition function Zq can be exactly calculated as 

ZQ = e-''^^lQ{yxt), (A2) 

where Iq is the modified Bessel function. If the form of G{9) is given, Gq can be also 
calculated exactly. For example, let us consider the case that G{9) = {9\Q\9) or {9\Q'^\9). 
From {9\Q\9) = -iVdE/d9 = -iVxtsm9, we have 

{Qh = ^ £|^X.sin^e--(^)--« = Vx/-^^^^^^^^^^. (A3) 

Using a formula for the Bessel function In-i{z) — /„+i(z) = {2n/ z)In{z), the correct result 
{Q)q = Q is reproduced. Similarly, using 

(^l^'l^) = + ^^1^1^)' = VXtCos9 - iVxtsm9)', (A4) 

we obtain an expression 
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which leads to the correct result {Q'^)q = Q^. 

We now consider a general G{9). We will calculate Gq, using the expansion of G{9) 

Gg = G(0) + ^-pr), (A6) 

n=l 



where 



This expansion can be evaluated using the connected part 



log 



and the formulas 



(A7) 



(X^) = (X='), + 3(X2),(X) + (X)^ 

(X^) = (X^), + 4(X3),(X) + 3((X2),)2 + 6(X2),(X)2 + (X)^ 



(A9) 
(AlO) 
(All) 
(A12) 



with X = ^ in our case. We evaluate the Q derivatives in the large V limit, using the 
asymptotic expansion of the Bessel function 



oo 



n=0 



(2.) 



where z = XtV and 



^ (4Q^-l^)(4g^-3^)---(4g^-(2n-in 

(Q, n) = , (g, 0) = 1. 

Using these the partition function is written as 

1 



V2 



TTZ 



i_ + + 



2z {2zy 



and we obtain ignoring 0{V ^) contributions 



log Zq = - log VQjtz - — 



(g,i) 



{QAf (g,2) 

2(2^) 2z 



logv^-l(g^-i) (i + l 



(A13) 



(A14) 



(A15) 



(A16) 
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from which 



logZQ = + (A17) 



dQ " ^ z \ 2z 
„ 1 1 



Therefore, 



^^^.ogZ« = --^l + -j. (A18) 
^logZg = 0, n = 3,4,5,--- . (A19) 



{0)q = 4 + ) ' (A20) 



= ^^ + 0{z~'), (A22) 
(^')q = 4 + (A23) 



Collecting these the final result becomes 



Gq = G(0) + G(^nO)^(l + ^ - + ^^'^^^^S^ ^^'^^^ 



+G«(0)z^(l + + G^^'^W^TT^ (A25) 



Noticing that C4 = — Xt for the present case, this result completely agrees with the previous 
result (!27l) within the 0(y~^) errors. This demonstrates that the saddle point expansion 
reproduces correct results, up to exponentially small corrections. 
The leading order of the chiral perturbation theory gives 

E{9) = N]xt [1 - cos(^/A^/)] . (A26) 

In this case, we can extend the integration range to —Nj-k < 9 < Nj-k by ignoring exponen- 
tially small contribution. After performing 9 integration exactly, we obtain 

Zq = e-'I^iz) + 0(e-"[^-™'^("/^/)l) (A27) 

where z = VNjXt and Q = NjQ. Noticing that C4 = —Xt/Nj in this case, it is easy to 
see that the saddle point expansion agrees with exact results order by order in the 1/V 
expansion, as long as exponentially small corrections are ignored. 
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APPENDIX B: CONTRIBUTIONS AT V'^ 



1. CP-even contributions 



For 9c = 0{V CP-even contributions become 



G 



Q 



2x 



3x^ 



4x3 



Using the large separation limit G = G^^) = g, G^^^ = 2Axt and G^^) = 480xtC4 for G = 
{ijj{xi)uj{x2)uj{xz)uj{x4)) , we obtain 



(u;(xi)u;(x2)u;(x3)t^(x4)) 



Q 



1^2 



- Q\t) 



o{v- 



(B2) 



2. CP-odd contributions 



For 9c = 0{y CP-odd contributions become 



^odd 



39r 



3! 2x \ 



2x2 y 4!x2 
5^(3) y 



12x3 



2x2 V 3x2 y g^;' 

Using the large separation limit G^^-* = 0, G^^) = q(5) _ ^gg^tQ for G 

(co'(xi)cj(x2)Lt;(x3)), we obtain 

7c4 



(B3) 



(u;(xi)u;(x2)c^(x3))( 



\/2 



3Xt + 



(B4) 
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